Suppression of Kondo screening by the Dicke effect in multiple quantum dots 
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The interplay between the coupling of an interacting quantum dot to a conduction band and its 
connection to localized levels has been studied in a triple quantum dot arrangement. The electronic 
Dicke effect, resulting from quasi-resonant states of two side-coupled non-interacting quantum dots, 
is found to produce important effects on the Kondo resonance of the interacting dot. We study 
in detail the Kondo regime of the system by applying a numerical renormalization group analysis 
to a finite-t/ multi-impurity Anderson Hamiltonian model. We find an extreme narrowing of the 
Kondo resonance, as the single-particle levels of the side dots are tuned towards the Fermi level and 
"squeeze" the Kondo resonance, accompanied by a strong drop in the Kondo temperature, due to 
the presence of a supertunneling state. Further, we show that the Kondo temperature vanishes in 
the limit of the Dicke effect of the structure. By analyzing the magnetic moment and entropy of 
the three-dot cluster versus temperature, we identify a different local singlet that competes with the 
Kondo state, resulting in the eventual suppression of the Kondo temperature and strongly affecting 
the spin correlations of the structure. We further show that system asymmetries in couplings, level 
structure or due to Coulomb interactions, result in interesting changes in the spectral function near 
the Fermi level. These strongly affect the Kondo temperature and the linear conductance of the 
system. 

PACS numbers: 73.63. Kv, 72.15.Qm, 72.10.Fk, 75.75.-c 



I. INTRODUCTION 

Quantum dots (QDs) have played a prominent role in 
the investigation of the Kondo problem physics in recent 
yearsji~— as they allow a systematic and well-controlled 
variation of structure parameters, with the consequent 
exploration of electronic correlations in this paramount 
many-body state. This exploration is of course possible 
by being able to experimentally tune the relevant param- 
eters in the system over rather wide ranges. Moreover, 
in a feature that is unique to these systems, QDs also 
facilitate the incorporation and study of quantum coher- 
ence effects in the Kondo state, and its interplay with 
structure resonances and field- induced phase shifts. In 
fact, the Kondo effect has been studied in multi-quantum 
dot systems, including double^— and triple dotsr 3 - - — in 
different geometrical arrangements, motivated in part by 
the richness of the Kondo physics in the presence of lo- 
calized levels, as well as by their potential application as 
spin filters. 

One of the main signatures of the Kondo state is the 
enhancement of the quasi-particle density of states at the 
Fermi level; this Kondo resonance is accessible in elec- 
tronic transport experiments, as it opens an additional 
transport channel, which is readily seen in the differ- 
ential conductance of these structures (typically in the 
zero-bias limit). In the case of multiple QD geometries, 
the Kondo and other single-particle-like resonances pro- 
vide different transport channels which can even interfere 
with one another. As a result, these structures offer the 
interesting possibility of controlling the transport proper- 



ties by exploiting quantum interference in the electronic 
propagation, allowing the study of scattering phenom- 
ena such as the well-known FancS and Aharonov-Bohm 
effects ; 23 ' 24 competing with the Kondo effect. For ex- 
ample, the interference of Kondo and Fano resonances 
in suitably designed structures has been shown recently 
in beautiful experiments, giving rise to complex con- 
ductance features! 11 ' 25 These results demonstrate that 
consideration of multiple coherent scattering of traveling 
electronic waves, both through single-particle as well as 
many-body resonances, is crucial for the understanding 
of the resulting conductance features. In this regard, be- 
ing able to tune the structure parameters of a multiple 
QD system, as well as its Kondo state, provides a unique 
arena to study correlations and coherence effects in a 
controllable manner. 

It is clear that the features of the Kondo resonance 
near the Fermi level significantly affect the conductance 
of the system, sometimes in quite subtle ways. For ex- 
ample, the shape of the density of states in the leads 
near the FL has been shown to produce strong modi- 
fications of the Kondo resonance and characteristic en- 
ergy (the Kondo temperature). These modifications can 
arise from intrinsic global properties ; 26 ' 27 or from the lo- 
cal environment in the vicinity of the active QD^ We 
are here especially interested in resonances and modifica- 
tions due to the electronic version of the Dicke effect^ 
The latter is the electronic analogue of the well-known 
Dicke effect in quantum optics, which takes place in 
the spontaneous emission of closely-linked atoms lying in 
the same environment (within one characteristic wave- 
length of each other) 2& In the electronic case, the de- 
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cay rates (level broadenings) are produced by the cou- 
plings between localized levels and a conduction chan- 
nel, and their close proximity and effective coupling gives 
rise to effectively fast (super-tunneling) and slow (sub- 
tunneling) modes i^2r— Interestingly, this coherent single- 
particle physics results in strong changes of the Kondo 
screening, as we will discuss below, once Coulomb inter- 
actions are fully taken into account. 

Although other multidot geometries would exhibit sim- 
ilar physics, a specific configuration to study this effect 
consists of three QDs where two large (essentially non- 
interacting) dots are attached laterally to a central dot 
that is embedded between current leads, as schematically 
shown in Fig. [TJ This cross-bar structure is reminiscent 
of quantum wave guides with a resonator cavity— In our 
case, the side QDs act as scattering centers to the main 
transport channel, and compete with the Kondo effect of 
the central doti^ Similar configurations have been stud- 
ied before, and recent work by Trocha and Barnas using 
a slave boson mean field approach has identified interest- 
ing regimes^ They studied the interplay between Kondo 
and Dicke resonances and considered the limit of infinite 
Coulomb repulsion (U — > oo), which suppresses virtual 
processes involving the doubly occupied state in the inter- 
acting dot. Additionally, the mean-field approximation 
adopted in the auxiliary boson fields neglects correlations 
due to charge fluctuations involving the empty state of 
the dot. All of these processes may be expected to be 
important in providing a full description of the Kondo 
resonance, especially in competition with the Dicke ef- 
fect here, as we indeed show in this work. Our numerical 
renormalization group (NRG) approach;^— provides a 
reliable description that incorporates all charge and spin 
fluctuations in the problem, it allows the study of finite U 
values-corresponding to real QD experiments-and gives 
us a systematic way to study the subtle interference ef- 
fects inherent in this geometry. 

In this work, we explore the Kondo state that appears 
when the levels of the non-interacting dots are brought 
together both symmetrically and asymmetrically about 
the Fermi level. In both situations, the Kondo reso- 
nance narrows drastically and is "squeezed" between the 
two single-particle levels, suggesting a drop in the Kondo 
temperature of the system, as described recently in the 
symmetric case^ We show that indeed the Kondo tem- 
perature drops, as the onset of the Dicke effect involving 
the interacting dot totally suppresses the Kondo state 
and results in vanishing Kondo temperature. This be- 
havior, due to the enhancement of spin-spin correlations 
in the QDs orbitals, results in the formation of a local 
singlet, which decouples from the current leads. As the 
system moves away from a fully degenerate Dicke con- 
figuration, the Kondo temperature is non-zero but has 
a non-monotonic dependence on structure parameters, 
behavior which is only brought out by our numerical 
renormalization group approach. This optimization of 
the Kondo effect in a structure may have interesting ap- 
plications in a real experimental system. 




FIG. 1. (color online) Schematic representation of the three- 
dot system. QD1 is considered to be strongly interacting (U 7^ 
0), whereas QD2 and QD3 are assumed to be effectively non- 
interacting (£/' small or vanishing-see text). 



We will further show that the competition of Kondo 
and Dicke effects modifies excited states near the Kondo 
regime as well, resulting in unexpected changes in the 
thermodynamic and transport properties of the system. 
Interestingly, the local (non-Kondo) singlet state men- 
tioned above emerges due to a strong coupling between 
the interacting QD and the super-tunneling state. A 
crossover between these two configurations can be tuned 
by varying the coupling between the QDs. 

We also describe the effects of asymmetries in the sys- 
tem, including different intcrdot couplings, as well as en- 
ergy levels and finite Coulomb interactions in all the dots 
of the structure. We show that these structural asymme- 
tries introduce interesting particle-hole asymmetries in 
the problem, resulting in strong changes in the Kondo 
temperatures and the linear conductance of the system. 



II. MODEL AND THEORETICAL APPROACH 

The system consists of an interacting QD (labeled 
QD1) coupled to two current leads (L and R) and to 
two effectively non-interacting QDs (QD2 and QD3) (see 
Fig. [1]); one could also think of QD2 and QD3 as being 
close to a Coulomb blockade peak, so that they conduct 
and their behavior is single-resonance like, as charging 
effects are typically small when conducting. The three- 
dot structure is described by a multi-impurity Anderson 
Hamiltonian, H = H dots + Hi eads + H T , with 
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where c\ a {c\ a ) is the operator that creates (annihilates) 
an electron with energy £j (i = 1,2,3) and spin a in 
the respective QD and c\ ka is the corresponding fermion 
operator for the leads £ — R, L, with energy e^, where k 
is the momentum quantum number of the free conduction 
electrons. The second term in Eq. (JTJ) accounts for the 
Coulomb repulsion in the doubly occupied state of the 
dots. The hopping amplitudes V, V12 and V13 couple 
the interacting QD1 to the leads and to QD2 and QD3, 
respectively. For simplicity, unless stated otherwise (see 
Sec. HUH), we will take V 12 = V 13 = V, and U 1 = U^0, 
while U 2 = U 3 = 0. 

The characteristic properties of the Kondo physics of 
the system will be probed through thermodynamic quan- 
tities, namely the local magnetic moment and its contri- 
bution to the entropy, which can be readily obtained from 
the NRG procedure. To calculate dynamical quantities, 
such as the spectral functions and the conductance of the 
system, one needs to calculate the local retarded Green's 
function (GF) at the interacting dot, which is defined in 
the standard form, 

GnM = ((ci .;4a»o> 

e- <UT «c 1(T (t);4(f)»dr J (4) 

where r = t — t' , and 

«c lff (t);4(tO)) = -ie(i-f)<[ci«r(*),cL(* , )]+) (5) 

is the double-time Green's function, with 6(r) the usual 
step function, and [■ • • ,•••]+ the anticommutator. Once 
Gl^w) is known, the spectral function p\ a {u)) [or density 
of states (DOS)] of the interacting dot can be obtained 
from the relation 
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A. Non-interacting analysis 

In the non-interacting limit (Ui — > 0) the GF G^ cr (w) 
can be easily obtained by the equation-of-motion method, 
leading to an exact closed expression, 
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In these expressions, the frequency lo must be understood 
in its analytical continuation sense, lo — > lo + in, where 
77 is an infinitesimal. Notice that the poles at £2 and £3 
appearing in the self-energy describe localized states in 
the effective conduction band due to the presence of QD2 
and QD3. In the limit of V —> 0, the self-energy reduces 
to the single impurity case, 
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where 



where Ao = 2irV 2 /2D is the hybridization function for a 
flat density of states, and D is the half-bandwidth of the 
conduction electrons. Notice that the effect of the leads 
and the two side-coupled QDs is fully taken into account 
in Eq. <[7|) by the self-energy Ei(w), which provides a 
"structured" effective conduction band. The electrons 
localized in QD1 are then coupled via the effective density 
of states in the leads represented by Ai in Eq. ©. 

The spectral function, p\ = p^ = pi in this non- 
interacting case is depicted in Fig.[2]for t\ = 0, V' = 0.05, 
V — 0.1 and different values of S. We observe that by in- 
creasing S the width of the central peak increases, while 
the satellite peaks are just slightly shifted away from the 
Fermi level. The shift of the satellite peak position can 
be estimated by analyzing the poles of G^. Assuming 
the satellite peaks to be Lorentzian-shaped, we find that 
their positions are given by ±y/2V' 2 +6 2 . In the inset we 
show the peak positions as function of S. The continu- 
ous (black) curve corresponds to the analytical expression 
1.01\/2V /2 + S 2 , while symbols correspond to the values 
obtained from the curves in the main panel. The need for 
the factor 1.01 results from the assumption that the peak 
has a Lorentzian shape in the analytical estimate. On the 
other hand, close to zero energy (lo ~ 0) and for small 
detuning (S <C A ), the DOS can be written as a su- 
perposition of a symmetric Fano lincshape-arising from 
the super-tunneling Dickc state-of width T + cx Ao, and 
a Lorentzian line-arising from the subtunneling state- of 
width T_ tx (5 2 /Ao, respectively— This complex behav- 
ior of the DOS is due to the hybridization of the lateral 
dots through the continuum of the central dot and leads. 
This phenomenon is in close analogy to the Dicke ef- 
fect in quantum opticsi 2 ^ Here, the super-tunneling and 
subtunneling modes give rise a broad (r+) antircsonance 
and a sharp (r_) resonance, respectively. As we will 
discuss later, this effect also modifies dramatically the 
many-body Kondo state. 



Ei(w) = Ai(w) + iAi(w) (8) 

is the self-energy in the non-interacting case, which takes 
into account the effects of the QD1 contacts with the 
leads, and with QD2 and QD3. The self-energy can be 
written as 

Sl H = 2y^-i- + ^ 2 "-g~ £3 (9) 

^uj-e k ((J - e 2 )(Lo - e 3 ) 



B. Interacting central dot 

We will now consider the interesting case of interac- 
tions in the central dot only, so that Ui = U is finite, 
while U2 = U 3 = U' = 0. In this case, the system can 
also be mapped into a single QD coupled to an effec- 
tive conduction band. In this case the dressed GF of the 
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the problem considering all the dots composing the "im- 
purity region," so that the three-dot cluster is coupled to 
a flat conduction band. This requires a simple generaliza- 
tion of the NRG procedure: start with a cluster impurity 
described by a set of 64 basis states, on which one writes 
and diagonalizes the initial Hamiltonian, Hdots , while the 
current lead is characterized by a flat density of states. 
For these NRG "cluster impurity" runs, we typically keep 
1600 states and use a discretization parameter A = 3.5, 
which provides sufficient accuracy for the spectral func- 
tion in the scales of interest. 



III. THE ROLE OF ELECTRONIC 
INTERACTIONS 



FIG. 2. (color online) Spectral function vs. energy for V — 
0.1, V = 0.05 and ei = 0, for 8 = 0.001 (solid, black), 
8 — 0.01 (dashed, red), 8 = 0.02 (dot-dash, green), and 
8 — 0.03 (dot-dot-dash, blue). The inset shows the posi- 
tion of the right-side satellite peak as function of 8. Solid line 
corresponds to the analytical expression (see text) while (red) 
circles show peak positions taken from curves in main panel. 



interacting dot can be formally written as: 
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where is the proper self-energy which takes into 

account the influence of the lateral QDs, the current 
leads, as well as the Coulomb repulsion U. Unlike the 
non-interacting case, a closed expression for Ei'(w) is not 
possible, due to the many-body term in the Hamiltonian. 
However, G xi (lj) can be reliably obtained using the NRG 
procedure . 38 ' 39 For convenience, we write the GF in the 
Lchmann representation^ 
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where Z is the partition function in the canonical en- 
semble, |n) is an eigenstate of the Hamiltonian H, with 
corresponding eigenenergy E n , and (3 = (fc^T) -1 . The 
imaginary part of G^w) is calculated at zero tempera- 
ture (T — 0) directly from the NRG spectrum^ while 
the real part can be obtained from a Kramcrs-Kronig 
transformation^ 



C. Cluster impurity approach 



A. Particle-hole symmetric case 

Hereafter, we set D = 1 (typically the largest energy 
scale of the problem) as our energy unit. We further con- 
sider the Coulomb repulsion U = 0.5, and E\ = —U/2, 
which corresponds to a particle-hole symmetric system 
with interactions in QD1. In order to study the Dicke 
effect in this system we first focus on the particle-hole 
symmetric point of the entire system. To this end, we 
tune the on-site energies £2 and £3 of the non-interacting 
QDs (U 1 = 0) symmetrically displaced from the Fermi 
level, for which we set £3 = —£2 = S. With this choice 
of parameters, the last term of Eq. ([§]) shows two poles 
located at uj = ±5, which produce strong enhancement of 
the effective coupling Ai(w). This results in a dramatic 
distortion of the local interacting DOS projected in QD1, 
Pi(ui), with important implications on the properties of 
the system. We first fix 8 (= 0.01) and vary the intcrdot 
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As we will see in the next section, it is also useful to 
analyze the details of the states involving the three QDs, 
considering the system as a whole. This approach also 
proves essential when dealing with Coulomb interactions 
in all three QDs, as the approach of an effective density 
of states described in the previous section, is no longer 
applicable in the more general case. To that end, we solve 



FIG. 3. (color online) Spectral function vs. energy for V = 
0.1, 8 = 0.01 and coupling V = 0.0 (solid, black), V = 0.01 
(dashed, red), V = 0.025 (dotted, brown) and V = 0.04 
(short-dash, blue). The inset shows the same curves over a 
wider energy range. Here, U = 0.5, U' = 0, and D = 1. 

coupling V' . Figure [3] shows pi(w) (= pit = pii in Eq. 
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[6]), for several values of V , The solid (black) curve shows 
the case V = 0, which corresponds to the typical case of 
a QD coupled to a conduction band with a flat density of 
states. Notice the Kondo resonance at the Fermi level in 
addition to the two Hubbard peaks at ps ±U/2 (see inset; 
which give rise to Coulomb blockade peaks in the conduc- 
tance at appropriate gate voltages). The width of the 
Kondo peak is proportional to the Kondo temperature, 
which in this case isi£ T K = T K {V = 0) = 5.6 x 10~ 5 . 
As V increases, the low-energy structure of p\{w) is 
strongly modified. Two satellite peaks emerge symmet- 
rically about the central peak at to = 0. These peaks be- 
come more pronounced for larger V and the distance be- 
tween them increases somewhat more than 2V2V' 2 + S 2 , 
as one would expect from the non-interacting regime. 
The widths of the satellite peaks increase as V increases, 
similarly to the non-interacting case, where their width 
is given by ~ A0V' 2 . Notice also Fano antiresonances 
in p\ between the central and satellite peaks, located at 
bj « ±<5. These antiresonances result from destructive 
interference between the Kondo state of QD1 and the lo- 
calized levels in the lateral QDs<2i We will discuss this 
point in greater detail below. 

It is interesting to notice that the central Kondo res- 
onance peak width shows a non-monotonic dependence 
on V, first increasing with V , but decreasing after a 
particular value. This behavior, first identified in Ref.l35l 
from the slave boson approximation, suggests that the 
Kondo temperature has a similar behavior, although this 
suspicion was not explicitly verified. Using NRG, we find 
indeed that to be the case. In Fig. |U(a) we show Tk as 
function of V for 6 = 0.01, and obtain a clearly convex 
curve with a maximum value (~ 10T$-) at V ~ 0.04, 
in agreement with the behavior seen in the Kondo reso- 
nance in Fig. |21 We should emphasize that Tk for V' > 
is larger than for the isolated QD1, T K (see value above). 

The zero-energy resonance also depends strongly on 5. 
Figure [5] shows pi(u>) for V = 0.05 and several values 
of S. Notice that all curves show well-defined peaks and 
their positions remain nearly constant at either cj w 
or ui w ±\/2V' 2 + S 2 . However, the dips between the 
Kondo and the satellite peaks are less pronounced for 
larger 6, and their position shifts with 6, as can also be 
seen in the inset. The displacement of the dips toward 
the Fermi level as 5 decreases produces a "squeezing" of 
the Kondo peak at the Fermi energy. The most evident 
situation shown in Fig. [5] is the case of S — 0.001, when 
the Kondo resonance takes the form of a narrow and 
sharp peak. The behavior of Tk with 6 is shown in Fig. 
2Kb), and it is clearly non-monotonic as well. Notice that 
for V' = 0.05 we find a maximum value of Tk — 20T K 
at S ss 0.025. More importantly, the "squeezing" of the 
Kondo resonance as <5 — > is accompanied by a strong 
suppression of Tk, which vanishes completely as S — > 0. 
The above behavior is a clear signal of the influence of 
the subtunneling mode on the Kondo state. Indeed, one 
can understand these results in terms of the Dicke state 
physics. For small values of S, the subtunneling Dicke 
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FIG. 4. (color online) (a) Kondo temperature as function of 

l K — 



V for 8 = 0.01. Notice T K at V = is non-zero, Tg 



5.6 x 10 here, (b) Kondo temperature as function of S 
for V' — 0.05. Other parameters are as in Fig. [3] Notice 
Tk — > for S — > 0. The step structure seen here results from 
the NRG discretization of the energy, and corresponds also 
to the discrete temperature values for which we calculate the 
thermodynamics (see below). Continuous lines show smooth 
fit of NRG results as guide to the eye. 
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FIG. 5. (color online) Spectral function vs. energy for V = 
0.1, V = 0.05 and 6 = 0.04 (short-dash, blue), 8 = 0.03 
(dotted, brown), 8 = 0.02 (dashed, red), 8 — 0.01 (dot-dash, 
magenta) and 8 — 0.001 (solid, black). The inset shows a 
zoom of the region around u = 0. 



state dominates and the Kondo temperature decreases. 
For sufficiently large values of 5 the supertunneling Dicke 
state dominates and the Kondo temperature reaches its 
maximum value. For larger 5 values, however, the sys- 
tem goes out of the Dicke regime and the Kondo tem- 
perature drops^ We should also comment that the total 
suppression of the Kondo screening for 6 = is essen- 
tially a quantum phase transition, which results in the 
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QD cluster isolating itself from the current leads. Per- 
haps one way to view this suppression of Kondo screening 
is to think of a Kondo resonance in the central dot that 
interferes destructively with the single-particle-like reso- 
nances of the side-connected dots (themselves hybridized 
with the leads via the central dot). The fact that all these 
resonances appear at the same energy, lj = 0, results in 
the eventual suppression of the Kondo state for the entire 
cluster. Of course, this is only a qualitative argument, 
ultimately validated by the NRG results. The underly- 
ing physical mechanism for this destructive interference 
and resulting ground state is fascinating, as we proceed 
to describe. 

To provide a different perspective on the suppression 
of the Kondo screening (and vanishing Kondo tempera- 
ture), we study the thermodynamical properties of the 
system, namely the contribution of the impurity to the 
magnetic moment, /i 2 , and the entropy, S, from which 
one extracts information about the nature of the states 
the system visits when the parameters are varied. As 
usual, M 2 (T) = tf ot (T) - M 2 (T) and S(T) = S tot (T) - 
S (T), where X tot (T) and X (T) (with X = n 2 ,S) 
are thermodynamical quantities of the system calculated 
with and without impurity, respectively. The magnetic 
moment and the entropy can be written in terms of 
the magnetic susceptibility and the partition function as 
M 2 (T) = T X (T), and S = \ogZ (we set k B = 1). In order 
to study the effects of all three dots, we employ here the 
"cluster impurity" approach described in Sec. Ill CI 

Figures M, a ) and[7Ja) show, respectively, fi 2 and S as 
function of the NRG iteration energy (temperature T) 
for several values of V' and S. Before engaging in the 
discussion of the more complicated situations, lets us ex- 
amine the result for V = S = (small open circles, 
purple) , corresponding to a single impurity case, as QD2 
and QD3 have their level at the Fermi energy (6 = 0) 
but are decoupled from the rest of the system (V = 0) 
we include this curve as a reference for the more general 
cases. For T ^> U, all the curves show convergence to 
/i 2 = 3/8 and S = log(64), resulting from the free or- 
bital (FO) fixed point of the system. As T -> 0, DQ1 
undergoes a crossover into a Kondo singlet state and its 
magnetic moment is screened by the conduction band, 
while QD2 and QD3 remain in their free orbital states 
(since V = and (5 = 0), resulting in a net unscreened 
magnetic moment of 1/4 and entropy S = log(16) for the 
three-dot cluster. 

For V',5 7^ 0, different regimes can be found. In the 
low temperature limit (T < Tk) the system always flows 
into a Kondo state, where QD1 is in a Kondo singlet con- 
figuration with the reservoirs, while QD2 and QD3 are 
combined into a symmetric/antisymmetric pair of states, 
one of which is doubly occupied (and the other empty) 
at T = 0, as schematically shown in Fig. |U[a). This situ- 
ation is characterized by a magnetic moment completely 
screened by the conduction band (/x 2 = 0) and with en- 
tropy 5* = 0, as shown by all curves in Fig. [S] and [7] 
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FIG. 6. (color online) (a) Magnetic moment for the three-dot 
cluster as function of temperature for S = 0.01 and various 
values of V' . Open small circles (purple) show results for 
V' = S = 0. (b) Magnetic moment as function of temperature 
for V' — 0.075 and various values of 8. 
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FIG. 7. (color online) (a) Entropy of three-dot cluster as 
function of temperature for 5 = 0.01 and various values of V' . 
Open small circles (purple) show results for V' = 5 = 0. (b) 
Entropy of cluster as function of temperature for V — 0.075 
and various values of S. 



Other interesting situations occur when < T < 8^ 
For V' < 6, the physics can still be understood in terms 
of the local orbitals, in which case (since T < 5) QD2 
and QD3 are still locked in one doubly-occupied and 
one empty state. In contrast, QD1 is in the magnetic 
moment regime, contributing alone to the magnetic mo- 
ment /i 2 = 1/8 and to the entropy S = log(2). When 
V' > S, however, the system goes into a more intricate 
configuration of local molecular orbitals. This is better 
understood in terms of symmetric (super-tunneling) and 
anti-symmetric (sub-tunneling) combination of the QD2 



and QD3 local orbitals, defined as 



1 ( 

71 (C2 ° 
7i (C2 ° 



c 3o 



C30 



(13) 



where c +a and c_ 



annihilate an electron in the bond- 
ing and anti-bonding orbital, respectively, with energy 
e + = £- = (for e 3 = — e 2 = 8 as before). With this 
transformation, the Hamiltonian of the dots reads 

Hdots = ^ £i c \a c i<? + u n it n n 

a 

+V2V' [ c Lc +CT + h.c. 

a 

-5j2[claC-a + h.c]. (14) 
a 

The transformed system is schematically represented in 
Fig.HJb). Note that the orbital "+" and "— " arc coupled 
to each other by a matrix element —5 and only the orbital 
"+" couples to the QD1 via matrix element v2V'. On 
this transformed basis we can provide a more intuitive 
analysis of the thermodynamic quantities as follows. For 
T K <T < 5 and V ^> 6 the configuration of the system 
is represented in Fig.^c), where the QD1 orbital couples 
to the "+" orbital to form a "local singlet" (essentially 
decoupled from the leads), while the "— " orbital remains 
almost free (since 6 <C V). The contribution to the 
magnetic moment and entropy is then provided solely by 
the "— " orbital, resulting in // 2 = 1/8 and S = log(4), 
as clearly seen in the curves with (blue) square symbols 
in Fig. inja) and [3a). We should notice that the two 
curves for V' = (for S = and S = 0.01)-in each of the 
(a) panels of both figures-are identical for temperatures 
T > 5 = 0.01, since for high T the QD2 and QD3 levels 
at ±<5 are thermally accessible (which is no longer the 
case when T <C S). 

In Fig. IHb) and Gib) we fix V = 0.075 and plot the 
magnetic moment and entropy for different values of 5. 
For 6 = 0.02 <C V, the situation described above is 
favored, but when 5 —> V' , the plateaus in the suscepti- 
bility and entropy shift towards the curves with empty 
(black) circles of panel (a) in the respective figures. The 
system is then returning to the situation when the local 
orbitals are described independently, as shown in Fig. 
IHd). Notice the Kondo singlet is not yet formed at these 
high temperatures (T > Tk)- These results reveal the 
establishment of a local antifcrromagncticaly correlated 
state, involving orbitals £1 and "+", emerging at temper- 
atures above Tk- 

In order to provide more evidence of the local singlet 
between QD1 and the "+" orbital, Fig. |H] depicts the 
spin-spin correlation value (Si • S+) (top panel); the cor- 
relation with the "— " orbital, (Si • S_) is shown in the 
bottom panel, as function of temperature for various val- 
ues of V . Notice that (Si • S+) is always negative (an- 
tifcrromagnctic) and has increasingly negative values for 




FIG. 8. (color online) Schematic representation of the various 
configuration of the systems, (a) Kondo singlet for <5 7^ and 
T < Tk- (b) Symmetric ("+") and anti-symmetric orbital 
("— ") combination, (c) Antiferromagnetic singlet formation 
between ei and "+" orbital, (d) Magnetic moment configu- 
ration for S ^ and T > T K . 



larger V' . In contrast, the correlation between the spin 
in QD1 and in the "— " orbital is an order of magnitude 
weaker, and totally vanishes for large V . It is interesting 
to notice that although the local correlation with the "+" 
orbital is large, it does not fully exhaust the anticipated 
singlet correlation of —3/4, even for the larger V' values, 
where it seems to saturate at ~ —1/2. 




V'=0.010 
=0.020 
=0.025 
=0.050 
=0.075 
=0.100 



10" 6 10" 5 10" 4 10" 3 10" 2 10" 1 
T/D 



10" 



FIG. 9. (color online) Spin-spin correlations as function of T 
for S — 0.01 and various values of V' . Lower panel proves that 
the spin correlation of QD1 with "— " is vanishingly small. In 
contrast, upper panel exhibits strong antiferromagnetic cor- 
relation with the "+" orbital with increasing V' . 



B. Effects of particle- hole asymmetry 

In a typical experimental situation, asymmetry in the 
structure parameters is more likely to occur. It is there- 
fore of interest to analyze the stability and behavior of 
the results discussed above with respect to changes that 
more generally make the problem particle-hole asymmet- 
ric. We introduce asymmetry parameters k v and Kg that 
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FIG. 10. (color online) Spectral function vs. energy for asym- 
metric coupling to the side dots. Top panel (k„ = 0.025) and 
bottom panel (k s = 0.05) V = 0.1, V = 0.05 and 5 = 0.001 
(solid, black), 5 = 0.01 (dashed, red), S = 0.02 (dot-dash, 
green), and 8 = 0.03 (dot-dot-dash, blue). 



FIG. 11. (color online) Spectral function vs. energy for the 
asymmetric bare levels about the Fermi level, with asymme- 
tries as — 0.025 (top) and Kg = 0.050 (bottom). The other 
parameters are V = 0.1, 8 = 0.01 and V' = 0.01 (solid, black), 
V = 0.02 (dashed, red), and V' = 0.03 (dot-dash, green). 



allow for V12 7^ Viz, an d | £3 ] 7^ [ea | , respectively as: 

V u = V(1±k v ) (15) 

and 

e l = 5(l±K S ), (16) 

where the positive (negative) sign in the expressions 
above corresponds to the index 1 = 2 (3). The fully sym- 
metric situation previously described corresponds obvi- 
ously to the case k v = ks = 0. 

We first allow for asymmetry in the couplings to the 
side dots; Fig. llOl shows the results for weakly asymmetric 
couplings V12 and V13 for different values of S. The top 
(bottom) panel shows the case of K v = 0.025 (k v = 0.05). 
The asymmetry in the couplings results in a shift of the 
central peak away from the Fermi level, as well as its 
strong suppression (compare with Fig. [5]) . This effect is 
less pronounced for larger values of 8. Notice, for in- 
stance, that for 8 = 0.001 the central peak is almost 
completely suppressed. Also, comparing the curves in 
the top and bottom panels, we observe that for a given 
5 the shift and suppression of the peak is larger as n v 
increases. 

In Fig. [TT]we fix 8 = 0.01 and k v = (symmetric cou- 
plings to central dot), and show the spectral function as 
function of energy for two values of ks and different val- 
ues of V . In the top (bottom) panel we have Ks = 0.025 
(ks = 0.05). We notice that for small V (for instance 
V = 0.01), there is a strong distortion in the spectral 
function near the Fermi level for both values of ks . How- 
ever this distortion is much weaker when V' increases, 
so that for V' = 0.02 and 0.03, the asymmetry in the 
spectral function is nearly absent. 

One would expect, from seeing these results for the 
spectral functions, that the particle-hole asymmetry in 



the system would introduce non-monotonic changes in 
the Kondo temperature. In Fig. [12] we illustrate this 
point. Panel (b) shows the associated Kondo tempera- 
ture for the system parameters of Fig. [TUk . where the 
coupling to the side dots is asymmetric (with a value of 
k v = 0.025). As suggested by the increasing width of 
the zero-energy peak in Fig. [T0a . Tk increases with 5, 
reestablishing the Kondo state that was absent in the 
full Dicke regime (for 8 = 0), while larger values even- 
tually reduce Tk, although not to zero. Similarly, the 
asymmetry in side dot level location illustrated in Fig. 
II lb ,, results in a rapid increase of Tk with V' , shown in 
panel (a), so that the system has an order of magnitude 
larger Kondo temperature for V' ~ 0.05, before dropping 
for larger V' values. 

An additional important source of particle-hole asym- 
metry is the bare level position of QD1. In order to 
explore how it affects our results, we allow the level ei 
to be different from —U/2. For convenience, we keep 
ei = —0.25 and change U, which allows not only to ex- 
plore the effect of the particle-hole asymmetry but also 
the effect of larger values of U. Figure [T3l shows the spec- 
tral function vs. energy for 8 = 0.01 and various values of 
U . The curve for U = 0.5, corresponding to the symmet- 
ric case studied in Sec. MI A\ exhibits two peaks symmet- 
rically placed about the Fermi level, as well as Hubbard 
bands located near t\ and ei + U, as can be seen in the 
inset of the figure. As U increases, one clearly observes 
that the curves no longer display particlc-holc symmetry. 
We observe a decreasing height of the central peak, as 
well as a minor shift in position, while its width remains 
almost unchanged with increasing U. More interesting, 
we notice that the satellite peak below the Fermi level 
is suppressed and shifted toward oj = 0, while the corre- 
sponding peak on the positive side is nearly unchanged. 
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FIG. 12. (color online) Kondo temperature for asymmetric 
three-dot systems. Panel (a) shows Tk (in units of T K - see 
caption Fig. [4| for system of Fig. Illb . Panel (b) shows Tk 
for system of Fig. HOfa . Notice strong enhancement of Tk in 
both cases with increasing parameter. 
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FIG. 13. (color online) Spectral function vs. energy for V = 
0.1, V = 0.05, S = 0.01 and U = 0.5 (solid, black), U = 
1.0 (dashed, red), U = 1.6 (dot-dash, green), and U = 2.0 
(dot-dot-dash, blue) (solid, black). The inset shows the same 
curves for the entire region within the conduction band. 



This is a clear consequence of breaking particle- hole sym- 
metry in the problem. 

We have also investigated the role of non-vanishing 
Coulomb interactions in the side dots QD2 and QD3. 
To that effect, the problem is addressed using NRG for 
the three quantum dot complex treated as a cluster im- 
purity, as described in Sec. Ill CI and with the inclusion 
of a non-zero Coulomb interaction U% = U3 = U'. Fig- 
ure [14] presents the spectral function for different values 
of U', for QD2 and QD3 levels placed asymmetrically 
about the Fermi level (6 = 0.01). For U' — the system 
is particle- hole symmetric (as e\ = — Z7/2), as also shown 



FIG. 14. (color online) Spectral function for the three dot 
cluster with non-zero Coulomb interaction in all dots for en- 
ergies close to the Fermi level. The solid (black) curve with 
U' = represents the same system as in Fig. [3j S = 0.01, 
V' = 0.05, V = 0.1, and U — 0.5. As U' increases, the asym- 
metry is evident, strongly affecting the Kondo peak near zero 
energy. 



in Fig. [3] For increasing U', however, the asymmetry 
in the spectral function is evident, strongly shifting and 
suppressing the Kondo resonance near the Fermi level. 
This indicates the important role that Coulomb interac- 
tions in the nearby dots have in the overall correlations 
present in the system, and suggests that the observation 
of the Dicke effect is indeed a delicate undertaking. 



C. Conductance 

We now turn our attention to the effect on the con- 
ductance of the system of the various features dis- 
cussed above. It is clear that the details of the den- 
sity of states at the Fermi level are directly connected 
to the zero-temperature transport properties in linear 
response. In the zero-bias limit, the conductance of 
the system can be obtained by the Landauer formula, 
G = (2e 2 /h) -Im[T(w)] (df(ui)/du) dui, where f (ui) 
is the Fermi distribution function and, in the present 
case, T(ui) = AoGdd(w). At T = the conductance 
reduces to G = (2e 2 /h)TrAopi(Ep). From the density of 
states in Fig. [3] and O we observe that for all values of 
V and S > 0, the DOS always exhibits a Kondo peak at 
the Fermi level (which becomes very narrow and sharp 
as S — > 0). Since the conductance of the system is pro- 
portional to the density of states of DQ1, this results 
in unitary conductance (2e 2 /h) for all values of V' and 
8 > 0, discontinuously dropping to zero at S = 0. The 
above behavior is explained by the complete localization 
of the subtunneling state for 6 = 0. In this case the state 
has no projection on the central quantum dot and there- 
fore it does not contribute to the conduction. Similarly, 
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FIG. 15. (color online) Conductance for asymmetric systems 
corresponding to Fig. 1121 for each of the two panels there. 



one would expect vanishing conductance for small bias 
as long as 5 = 0.— 

An interesting result of the particle- hole asymmetry, as 
we have seen in the previous section, is to shift the Kondo 
resonance away from the Fermi level, as well as to change 
its amplitude. This has direct impact on the linear con- 
ductance of the structure. Figure [15] illustrates this be- 
havior for different asymmetries (notice in fact that al- 
though k v and n$ are constant in this figure, increasing 
V or 8 value increases the asymmetry of the structure). 
In both cases, the conductance is indeed found to change 
drastically with the appropriate parameter. It is interest- 
ing that asymmetry in the QD2 and QD3 level position 
(for Kg = 0.025) and increasing dot coupling V' results 
in a strong suppression of the linear conductance, even 
as the Kondo temperature of the system (Fig. [T2"a ) is 
increasing. On the other hand, the enhancement of Tx 
that we see for the coupling asymmetry (k v = 0.025) in 
Fig. 112b . does produce a rapid increase in conductance 
for larger 8 values. 



IV. CONCLUSIONS 

In summary, we have investigated the interplay of 
Dicke and Kondo effects in a central strongly interacting 
quantum dot coupled to a conduction band and two local- 
ized levels provided by additional nearly non-interacting 
dots. We study in detail the Kondo regime of the sys- 
tem, by applying a numerical renormalization group anal- 
ysis to a finite- U multi-impurity Anderson Hamiltonian 
model. We find that the system displays a "squeezed" 
Kondo resonance as we tune the single particle levels 
of the side-coupled dots toward the Fermi level. In the 
strict Dicke limit of degenerate single-particle levels, the 
Kondo singlet disappears, as its Kondo temperature van- 
ishes. This behavior is found to be closely related to the 
formation of a local singlet state involving the (coupled) 
orbitals of the interacting dot and the super-tunneling 
state (the symmetric combination of the single particle 
levels of the two side coupled dots). This local singlet 
results in the suppression of the Kondo screening, which 
is detrimental to the conductance of the system. We 
have further explored the consequences of system asym- 
metries, introduced by different couplings, level locations 
and Coulomb interactions. These are found to strongly 
affect the spectral density and the conductance through 
the structure. We believe that our study not only clarifies 
the relevant picture of an interacting system coupled to 
side quantum dots and conduction band, but it also pro- 
vides important guidelines for experimental realization of 
these ideas. 
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